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Abstract 



We examine charge/colour breaking along directions in supersymmetric field space which 
are F and D-flat. We catalogue the dangerous directions and include some new ones which 
have not previously been considered. Analytic expressions for the resulting constraints 
are provided which are valid for all patterns of supersymmetry breaking. As an example 
we consider a recently proposed pattern of supersymmetry breaking derived in Horava- 
Witten M-theory, and show that there is no choice of parameters for which the physical 
vacuum is a global minimum. 
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Supersymmetric models often have charge and/or colour breaking minima which can 
compete with the physical vacuum |], || . Excluding the regions of parameter space that 
lead to such minima results in bounds which can be severeQ. There are a number of 
reasons why supersymmetry suffers from this kind of metastability. There are many flat 
directions which are lifted only by our choice of supersymmetry breaking parameters. The 
higgs scalar which couples to the bottom quark, Hi, has the same quantum numbers as 
the leptons. Moreover, the mass-squared parameter of the higgs which couples to the top 
quark, H 2 , is negative at the weak scale (a requirement of electroweak symmetry breaking 
and one of the successes of supersymmetry). 

The resulting charge and colour breaking bounds can be roughly divided into three 
kinds. 

• D-Hat directions which develop a minimum due to large trilinear supersymmetry 
breaking terms. They are important at low scales and are usually referred to as 
CCB bounds. 

• F and D flat directions corresponding to a single gauge invariant. These are impor- 
tant when there are negative mass-squared terms at the GUT scale. One therefore 
needs to take account of the renormalisation group running, but they can be lifted 
by non-renormalisable terms in the superpotential. They are often referred to as 
unbounded from below (UFB) bounds ||. 

• F and D flat directions which correspond to a combination of gauge invariants 
involving H 2 . These can develop a minimum at intermediate scales due to the 
running H 2 mass even if all the mass-squareds are positive at the GUT scale. They 
cannot be lifted by non-renormalisable terms. They are also (rather confusingly) 
referred to as UFB bounds. 

The 3rd kind of bound will be of interest here since they can be very severe. For 
example, the Constrained MSSM (CMSSM) near its low tan/3 fixed point has almost 
half its parameter space excluded by them. Combined with dark matter constraints and 
experimental bounds, this is sufficient to exclude the model entirely 0. Unfortunately 
these bounds are also the most difficult to calculate because the renormalisation group 
running of the soft supersymmetry breaking parameters plays a complicated role in their 
determination. Specifically, the dangerous minimum forms where the running (negative) 
H 2 mass-squared is offset by the other (positive) mass-squared terms. 

Hence most work on UFB bounds is done numerically and the models considered 
typically have few parameters (such as the CMSSM which has only four). The reason 
for this is a practical one; to properly sample a parameter space one has to include, 
say, one low and one high representative of each variable. This becomes progressively 
more difficult as the dimensionality increases and for the most general MSSM (which has 
~ 100 parameters) it is clearly impossible. Even if it were possible to sample such a large 
parameter space, extracting any useful information would be extremely difficult. 

1 The alternative, to devise a cosmological mechanism to place us in the metastable physical vacuum, 
was discussed in Ref.[fJ, |j and will not be considered here. In this paper we will be concerned with the 
analytic determination of these bounds and shall therefore take a pragmatic approach; even if cosmology 
does place us in the physical vacuum we would, at the very least, like to be able to calculate if the present 
vacuum is unstable. 
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However, what does one do if a case of phenomenological interest (the real World for 
example) requires a less restricted set of parameters than that of the CMSSM? How do the 
UFB bounds restrict supersymmetry breaking in more general models? While determining 
this would very difficult numerically, it was shown in Ref.[|J that, by making a few simple 
approximations, we can make a reasonably accurate analytic approximation. In this 
letter we shall present the analytic expressions for the UFB bounds in the most general 
supersymmetry breaking scenarios. In particular we will show that there is essentially 
only one particular combination of supersymmetry breaking parameters which is subject 
to UFB bounds. This combination is 

2m\.. +ml-m 2 U33 -m 2 Qs3 , (1) 

where the notation for soft supersymmetry breaking parameters is as in Ref.[|J. 

We begin by ennumerating the dangerous directions, assuming the usual -R-parity 
invariant superpotential of the MSSM; 

Wmssm = huQH 2 U c + h D QHiD c + h E LE x E c + ^H X H 2 . (2) 

As we stated above, the dangerous F and D flat directions are those constructed from 
gauge invariants involving H 2 H § • The first example of this kind in the literature is (see 
Komatsu in Ref. [0) 

L { Q 3 D 3 ; H 2 U (3) 

where the suffices on matter superfields are generation indices. With the following choice 
of VEVs; 

h Q 2 = -a 2 n/h Dyi 
<^l 3 — dft 3 = a{i/ti£) 33 

Pi = a a/1 + a 2 /j,/h D33 , (4) 

the potential along this direction is F and D-flat, and depends only on the soft super- 
symmetry breaking terms; 

V = jl—a 2 (a 2 (ml + m 2 L J + m 2 Li . + m 2 d33 + m 2 Q J. (5) 

This is not quite, but is very close to, the deepest 'fully optimised' direction Two 
additional points. Note that in this case we can work in the basis in which the rele- 
vant Yukawa coupling is diagonal. Also, since the only severe bounds involve the third 
generation we shall restrict our attention to these cases. 

At large values of a ^> 1 the potential is governed by the first term. This can be 
negative at the weak scale because of the negative value of m\. The potential can therefore 
develop a charge and colour breaking minimum at a scale of few x /i/hoj and ensuring 
that this does not happen leads to the UFB boundQ. 

2 As discussed in Ref. , the traditional bound - that the physical vacuum is the global minimum - is 
very close to a sufficient condition - that the physical vacuum be the only minimum. The latter condition 
circumvents any cosmological questions concerning, for example, the tunneling rate between minima. 
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The coefficients in Eq.(|j) were chosen to cancel Fh x = dW/dH\. However there 
are other combinations of invariants which can do the same job; in fact L^H^ plus any 
invariant containing either L 3 E 3 or Q 3 D 3 can cancel Fh 1 - However the invariant should 
not contain Hi otherwise Fp, Fq and Fi would all be non-zero. The complete basis 
of gauge invariant monomials in Ref.|| then tells us that the following invariants are of 
possible relevance; 

LiL 3 E 3 i ^ 3 

LiQ 3 D 3 

Q^UjQ' 3 D 3 1 ± 3 

(Q'^XQ^XQ^s)^ +L perms 

(UMDjC&DjHz 5 3l + 6 3j + 

(U 1 U 2 D i )Q' 3 D j Q , l D n I ^ 3; 5 3i + <% + 5 m3 ^ 0, (6) 

where groups of three coloured fields are contracted with Levi-Cevita symbols. Note that 
in the above list, Q' is in the up-quark mass basis to ensure that Fjj 2 = 0. The first two 
directions have already been considered at some length in the literature. However this 
will be, to our knowledge, the first time that the other directions have been discussed. 

We can discard most of the new directions because they have non-zero Fq or Fjj. For 
example both the UUDQDH2 and UUDQDQD directions have non-zero Fq> 2 . This is 
enough to lift the radiative minima since the latter has a depth of ~ a^^^m^/h 2 where 

CLmin = O(10), but 

f q> 2 = h 2 c a 6 min fi 4 /h 4 b > a^ in ii 2 m 2 w /hl. (7) 

Clearly the only dangerous directions involve only the first generation up-quarks. For 
example, the Q' 1 UiL 3 E 3 direction is dangerous. Although it has non-zero F^, Fq> and 
Fh 2 , these F-terms are not enough to lift the radiatively induced minima since 

Similar consideration apply for the Q'UQ'D direction. Finally the Q'Q'Q'LLLE direc- 
tion has nothing apart from non-zero Fu 1 -terms if we choose VEVs along the hf^Vk an d 
U' Li D' L2 D' L , E L . i v i v j E 3 directions. 

The potentials for these operators (when combined with h^fk and neglecting the small 
Fy 1 , Fq,^ and Ffj 2 terms above) are 



LLE - -^-a 2 (a 2 (m 2 + m| fc J + v / ^TT(l-^)(m| ifc +mi fei ; 
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+ m\ u + m| 33 + m| 33 ; 



LQD -> -^-a 2 {a 2 {m 2 2 + m 2 L J + Va 2 + 1(1 - 5 ik ){m 2 Lik + m 2 L J 

n D 33 

+ + ml + m 2 , ■ 
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11 2 

r\UT TP 2/2/2, 2\, 2 1 2 1 2 , 2\ 

QULE -»■ -3— a (a (m 2 + m Lfe J + m Q , + m Uu + m Lg3 + m E J 
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QUQD -> ^a 2 (a 2 (ml + mlj +ml, ii +ml ii +ml, 3s +m 2 D33 ) 
QQQLLLE - » a 2 (a 2 (ml + m|J + (1 - <f y )(mL + m | ) 



+ Va 2 + l(l - <Sy)(mL + mi ) + Va 2 + 1(1 - 5 ik )(m 2 Lik + m 2 L J 



-m 



Q11 + m Q 22 + m L + m L + m L + m L + m L + m i 33 + m U) 



where the /i£) 33 should be the Yukawa coupling taken in the up-quark diagonal basis for 
the QUQD direction. At first glance the last direction seem to be safe because of the large 
number of mass-squared terms. However, it should be stressed that these parameters are 
not physical mass-squareds and they can therefore be negative (as is often the case in 
weakly coupled string models for example). As we shall see however, the renormalisation 
group running does tend to make them positive towards low scales so that they indeed 
tend to be less important. The off-diagonal bounds are likely to be unimportant when 
the relevant mass-squared parameters are positive, so henceforth we shall consider only 
the diagonal ones. 

In order to obtain the UFB bound we now need to take account of the renormalisation 
group running of the mass-squared parameters between the weak and GUT scales. An 
accurate analytic method was introduced in Ref.|4[] which we shall briefly recap. 

First we need to make some approximations and assumptions; they are 

• Neglect two loop effects in the renormalization group running of parameters (here 
we include hypercharge contributions). 

• Of the Yukawas, keep only h t = hu 33 in the runnning - neglect the bottom and tau 
Yukawas (i.e. tan/3 < 30) and mixing. 

• Assume that the gaugino masses are degenerate (M 3 = M 2 = Mi = M 1 / 2 ) at the 
scale where the gauge couplings unify (which we shall refer to as the GUT scale) 

These approximations allow us to obtain analytic solutions for the running parameters 
in closed form (which were presented in Ref.[Q). In Ref.|Q it was shown that the final 
analytic expressions for the UFB bounds should be accurate to ~ 15%. 

Next, we assume that the largest mass, and therefore the appropriate scale to evaluate 
the parameters at is = hjj.^h^). In the above potentials, (h®) = —a 2 fj,/hD 33 ,E 33 so that 
the potentials are of the form 

V = ^i($A + B/b) (10) 

""^33 

where A = m 2 .^) +m 2 L ..(<j)), B is the other combination of mass-squared parameters (also 
evaluated at <f>) which appears in the potentials above, 

4> = <p/M GUT (11) 

and 

= M GUT h D:i3 = 2.9 x 10 12 / M GUT \ / m b (4>) \ / 200GeV \ 
[<P) h U33 pi hvjl.l V2 x 10 16 GeVy l2.5GeV/ I ) 1 ' 
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for /i£> 33 or 
6(0) = 



M GUT h E33 2.0 x 10 



12 



h U33 n 



h Uaa /l.l 



M GUT \ ( m T (0) \ / 200GeV \ 
2 x 10 16 GeV/ V1.78GeV; V J 



(13) 



for tiE 33 - The traditional UFB bound is saturated by V = V = 0; the non-trivial solution 
is therefore also a solution to V = V' = where 



V = 4>A + B/b. 



(14) 



Hence, only the parameter b enters the bound. In fact the bound always becomes more 
restrictive with increase in b, since this decreases the positive contribution to the potential. 
(The increase with unification scale has already been noted in numerical work, e.g. Casas 
et al of Ref.||.) To solve for the UFB bound we first define the point P where 



A{<j> v ) = 0. 

We can then expand about <p p and solve V(4>d) = V'{ 



(15) 



0; 



eB 

~Fb 



(16) 



Likewise, the sufficient condition is given by, (V""(0 C ) = V'(0 C ) = 0); 

e 3 / 2 B 



2A'b 

,"3/2, 



(17) 



The solutions for the running parameters including b and A 1 were given in Ref . 0] in terms 
of the three parameters with quasi-fixed points; 



R 

A t 
3M 2 



°u 33 

m U 3 3 +m Q33 +m 2- 



Using these solutions we can find the value of 4> p by solving Eq.(jrBj) or Eq.(|T^) either 
numerically or iteratively. Finally we obtain the bound by solving Eq.flIB]) (which we 
haven't yet used); 



A(<f> P ) 
where 







\M 2 + (raj + ml u - \M-) 



(2)> 



(19) 
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(n) 



a? 



1 



(20) 



and the 0-subscript indicates values at the GUT scale. Using Eqs.( p!6|JT9| ) we are now able 
to estimate the UFB bounds in any particular case of interest and to get a reasonable 
estimate for the most general case. 
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General models at the low tan j3 fixed point 

First consider the low tan/3 fixed point (where hu 33 is very large at the GUT scale). 
As a rule of thumb the UFB constraints are most restrictive here. This is because the 
dangerous minimum is driven by the negative mass-squared of H 2 which is in turn driven 
by the top quark Yukawa coupling. A larger top Yukawa makes m\ negative closer to the 
GUT scale and the other mass-squared terms have to be larger to overcome it; and the 
largest possible top Yukawa is at the fixed point. 

Here we see why the bounds involve the particular combination appearing in Eq. ([!]). 



From Eq.(19) the relevant bound is 
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(m 2 2 + ml i --M 2 : 



> -M 2 - MU°6® + ±g&) + _^(2) } . (21) 
o 2 1/iV 9 s 2 2 198 



The M 2 parameter is determined at any particular scale since we are at the fixed point. 
The only influence the B/b term can have is logarithmic, via the determination of the 
scale, (fi p , at which the RHS should be evaluated. Consider the LLE, LH 2 directions with 
the central value of b. Using the solutions of Ref. M this scale is given by 



2ne ((m 2 Lii +m 2 L33 + m 2 E J\ -36i 2y 



22) 



3Ho V «3#(3M 2 + AfjJ - 2a| 
where m = rn/Mi/2, which can be solved iteratively. The final bound is 

{2m 2 Li% +m 2 2 - m 2 U:j3 - m 2 Q J q > f(B\ ) (23) 



where 



B\o = (ml u + m\ m + m% 33 

2 



f(x) w 1.43 - 0.16s + 0.02x 2 . (24) 

The last equation is a fit to the bounds which is accurate for x £0. We stress that 
this result is valid for the MSSM with any pattern of supersymmetry breaking. In the 
CMSSM, where the scalar masses are degenerate and equal to m at the GUT scale, the 
inequality becomes 

m 2 >f(3m ) (25) 

and we find the familiar result 

m >1.05M 1/2 (26) 

which is a slight overestimate of the numerical results of Ref.p| (by 5 — 10%). 

We can also examine the dependence on the b parameter for this case. For example, 
the central value of b had p = 200 GeV. If we take fi = 500 GeV (or b = 7.7 x 10 11 ) we 
find 

f(x) w 1.20 - 0.14x + 0.02x 2 (27) 
and for the CMSSM we then find 

m >0.95M 1/2 . (28) 
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From this we can conclude that, given the accuracy of our other approximations, there is 
little point in determining the value of // by minimising the potential (which is fortunate 
since this would have made our task prohibitively difficult). Additionally Mqvt is as 
important as p, in the determination of the bound and we can reasonably set /i = M1/2 
for an estimate of the UFB bound. 

Before leaving the fixed point, there is an additional fairly weak bound on B\ . When 
this parameter becomes too negative it is unable to lift the minimum for any GUT scale 
value of m 2 L .. + m 2 . The resulting bound is 



B\ 



[m, 



+ m\ 33 + m 2 E33 



> -Q.7M( 



/2" 



(29) 



The values of the functions, f(x), and the above bounds are given in the table for all 
the directions in Eq.(^). As can be seen the constraints corresponding to the LLE, LH2 in- 
variants are the strongest. This had been noted before in numerical work on the CMSSM, 
but here we see it is a completely general result. 



LiJ 2 +Operator 


/(*) 




B > 


LLE 


1.43 - 0.16a; + 0.02x 2 


2.94 - 0.2x + 0.02x 2 


-0.7 


QDE 


1.08 -0.05s 


2.55 - 0.06a; 


-3.5 


QULE 


1.03 -0.04s 


2.49 - 0.05a; 


-4.5 


QUQD 


0.94- 0.022 


2.38 - 0.03a; 


-7.2 


QQQLLLE 


0.96 - 0.02a; 


2.39 - 0.03a; 


-6.8 



General models away from the fixed point 

At the fixed point, the bounds do not depend on M 2 \ or ^4t/ 33 |o = A . As we move away 
from the fixed point (to higher tan/?) we find that the dependence on these parameters 
increases. We now include this dependence. The distance from the fixed point is most 
conveniently expressed in terms of the parameter 

R 

P = ^QFP ( 30 ) 

where R® FP is the fixed point value of R which is simply a function of scale. 

The bounds can be evaluated as above from Eq.(|l9~D keeping the dependence on A 
and M 2 \ . In order to determine the scale <p p we can approximate A'(4>) with its fixed 
point value which is independent of A and M 2 | ; the actual value is close to the fixed 
point value even for quite large values of tan/5 and in any case the scale appears only 
logarithmically in the determination of the UFB bound. The bounds are found to be 

{2fh 2 Lii +m\- m 2 U33 - m 2 Q33 )\ o > f(B\ )+{p p -l) (g(B\ Q ) + 3M 2 | - Pp (l - A ) 2 ) , (31) 

where p p is the value of p at the scale (f> p . To a very good approximation we find that the 
value of p p is given by || 

- = l + -L = l + 3.17(sin 2 /3 - sm 2 (3 QFP ), (32) 
p p 2R 
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where Ro is the value of R at the GUT scale and tan/3^ FP = 1.78 in this approximation 
(but it should be stressed that the bounds we are evaluating are not dependent on its 
absolute value but rather the value of p p ). In Eq.(|3~TD, f(x) is the bound at the fixed 
point as before, and g(x) can be approximated in a similar manner. For the LLE, LH 2 
direction we find 

g{x) ^ 2.94- 0.20x + 0.02x 2 . (33) 

The functions for the remaining directions are given in the table above. We can see that 
A « M1/2 results in the weakest bounds. In addition the bounds asymptote at large 
tan (3 to the values with 

= 1 + 3.1700^,9^ " °' 57 ' (M) 

To get some estimate of the errors involved in neglecting the two loop and threshold 
effects, when they are included in the running t&n j3® FP m 1.6, which would give 

p p « 0.53. (35) 
Application: UFB bounds in M-theory 

Now let us apply this method to supersymmetry breaking by bulk moduli superfields in 
Horava- Witten M-theory |7|] . The supersymmetry breaking terms were derived in Ref . || , 
for the case where the dilaton, S, plus a single modulus field, T, are responsible for super- 
symmetry breaking (see Refs.|| for more on this and other methods of supersymmetry 
breaking in M-theory). The supersymmetry breaking can be parameterized by defining 
the F-terms as 

F s = V3m 3/2 (S + S) sin 9 

F T = v/3m 3/2 (T + T) cos0, (36) 

where m 3 / 2 is the gravitino mass and we have (as in Ref.||) set the phases and the 
tree-level vacuum energy density to zero. The structure of the resulting supersymmetry 
breaking is substantially different from the weakly coupled string, and depends on 

y = a(T + T)/(S + S) (37) 

where a(T + T) comes from a correction to the gauge kinetic function; 

Mi/2 V3 / . y \ 

— — = sinH -= COS V 

m 3 /2 1 + 2/ \ V3 J 



m 2 



y 12— sin 0+12 — cos — sin cos 



m\ j2 3 + y\ y \ 3 + yJ \ 3 + y J " 3 + y 

A ° Vs((-l + ^)sme + V3(-l + -^-)co S e). (3J 



"T.3/2 VV 3 + VJ V 3 + 2/, 

Now from Eq. fl3"T| ) we can see that in order to avoid new charge a colour breaking minima 
we should satisfy the inequality 

F(y, 6) = m 2 - f(m 2 ) + (1 - p p ) (g(3ml) + 3m 2 - p p (l - A ) 2 ) > 0. (39) 
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The inequality is most easily satisfied away from the fixed point so we choose the asymp- 
totic low value of p p = 0.53. Fig.(l) shows the resulting contour in the (y, 6) plane. 
At first sight this would seem to be a good result since there are some regions where 
the inequality is satisifed. However, in and between the two regions our approximations 
break down because the mass-squared parameters are negative at the GUT scale and the 
potential never becomes positive - i.e. there is no p p . Between the contours the potential 
is really unbounded from below with no minimum. The UFB bound is actually closest to 
being satisfied at 9 « 0.4, however no reasonable choice of parameters (/i or tan f3® FP for 
example) can lift the charge and colour breaking minimum. 

Thus we conclude that in this model there is no choice of parameters for which the 
physical vacuum is the global minimum. This could have been guessed from the fact 
that (as stated in Ref|§) the scalar mass parameters are always less than the gaugino 
masses. However showing it numerically would have been difficult. We should also add 
that the minima are only a problem when i?-parity is conserved. i?-parity violating models 
with this pattern of supersymmetry breaking are safe from UFB constraints provided the 
i?-parity violation is strong enough and in the correct Yukawa couplings |4| . 
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Figure 1: The F(y, 6) = contour in the (y, 6) plane. 
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